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Abstract. The well-known astrophysical resonant-reaction-rate (RRR) equations for an isolated narrow res- 
onance induced by the charged particles have been reexamined. The validity of those 'look reliable' as- 
sumptions used in deriving the classical analytic equations has been checked, and found that these analytic 
equations only hold for certain circumstances. It shows the customary definition of 'narrow' is inappropri- 
ate or ambiguous in some sense, and it awakes us not to use those analytic equations without caution. As 
a suggestion, it's better to use the broad-resonance equation to calculate the RRR numerically even for a 
narrow resonance. The present conclusion may influence some work in which the classical narrow-resonant 
equations were used for calculating the RRRs, especially at low stellar temperatures for those previously 
defined 'narrow' resonances. 

PACS. 25.40.Ny Resonance reactions - 26. Nuclear astrophysics 

1 Introduction models as input. Thereinto the reaction rate, which is de- 
fined as how fast a reaction takes place, is one of the most 

Nuclear astrophysics addresses some of the most remark- cr iti C al ones. It reflects most of the critical stellar features 

able questions in nature: What are the origins of the el- ( such as time sca i ei energy production as well as nucle- 

ements making life on earth possible? How did the sun, osynthesis of the elements, etc). Nuclear astrophysics has 

the solar system, and the stars form and evolve? and etc.. been developed for more than 50 years, a great number of 

In order to answer these questions, a huge amount of nu- stellar rea ction rates have been calculated by utilizing the 

clear physics information is required in the astrophysical classical equations based on the experimental data. How- 
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sical reaction-rate equations derived in the historic books written in a form of, i.e., Equ. (4.53) in [3J, 
or references, e.0., in [HUE]. , , / 8 \ 1/2 1 / E R 

((TV) ' ' 



E R exp(-^jJ (r(E) N BW dE. 



7t fi J (fcr) 3 / 2J 

This work is focusing on the validity of resonant-reaction- 
rate (hereafter referred to as RRR) equations for an isolated And thcn thfi integration of the Breit-Wigner cross sec- 
narrow resonance induced by the charged particles. Actu- ^ yieldgj for ft narrQw rcsonancc with ( A_ )2 „ ( A.^ 
ally these equations have been derived in many books or (hereafter referred to as « the SEC OND assumption") and 
references. As an example, let's thumb through the famous w[th ncgligible energy dependenC e f the partial and to- 
Book " Cauldrons m the Cosmos" [3J. To make readers tal widthg (hereafter referred to as « the XHIRD assump- 
easy to follow, a similar deriving procedure and the same ti on >n j e jrjq U (4 54) m [3] 
notations used in [3J are being adopted here. 

It is well-known the stellar reaction rate per particle 
pair (av) can be written as (Equ. (4.1) in [3J) 

1/2 
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(E-E R y + {r/2) 

= 2tt 2 (A)^^. (5 ) 
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Where u>j=lu r "p b is defined as the resonance strength. 



• ) 1 - rn ) 3 / 2 Jo a ( E ' Eex v( kT^) dE ' Finally, the RRR per particle pair for an isolated narrow 
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resonance can be written as 



As for an isolated resonance the cross section can be ex- / \ 3/2 



2, / Er 



pressed by a famous Breit-Wigner formula as (Equ. (4.52) 



m jjj^ This classical analytic Eq. [6] is so famous that it has been 

utilized in a great number of nuclear astrophysics work. 

<j {E)bw = ^^2^' [2J\ + 1)(2J2 + 1) After a careful survey, we find the definition and bor- 

x - - (2) der of narrow and broad resonances in the previous liter- 

{E-E R f + {r/2f {1 

ature are not quite clear. For instance, a resonance with 
With the knowledge of the energy dependence of the cross r /E R >lQ % called a broad resonance [3J implies that a 
section for a narrow resonance (r«£ fl ), Equ. Q] of the resonance with r/E R <W % can be treated as a narrow 
stellar RRR per particle pair can be replaced by ig it appropriate? In the present workj we will exam _ 

/ Q \ 1/2 roo , j-, x ine the validity of all three assumptions mentioned above. 

<™>=y WW l ^)^ex P (--)^. 

(3) 

2 Validity of Assumptions 

Assuming the Maxwell-Boltzmann function Eexp(—E/kT) 

changing very little over the resonance region (hereafter The FIRST assumption Chronically a resonance for which 

referred to as "the FIRST assumption"), Equ. [3J can be r/E R <10 % was treated as a narrow resonance (e.g., see 
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Fig. 1. Amplitudes of relevant functions used in the equations. 



The SECOND assumption The de Broglie wavelength 
squared (t^:) 2 ( = ^g) is inversely proportional to the center- 
of-mass energy E, and therefore this assumption (^) 2 ~ 
means E~Er over the resonance region. It's a reasonable 
assumption for a resonance with r/En<\0 %. 

The THIRD assumption The partial width for a charged 
particle a (usually either proton or a-particle) is strongly 



The units are included in the square brackets. See text for energy dependent because the projectile must penetrate 



details. 



the Coulomb barrier, and can be defined as (Equ. (4.65) 
in ®) 



2Ti [2E\ 

Hn \ U 



[3]), here it shows this definition of 'narrow' is inappropri- 
ate especially at low temperatures. If one defines the varia- 
tion 5 M b of the Maxwell-Boltzmann function Eexp(-E/kT) where the penetration factor can be calculated by i^= (F ^ G 2 } 0] 
over the resonance region as 



. . E R + r/2 ( E R + r/2-E R 
Omb = f exp 



Er 



kT 



with the regular and irregular Coulomb wave functions 
Fi, Gg. The quantity 0j, called the dimensionless reduced 
width of a nuclear state, is generally determined experi- 
mentally and contains the nuclear structure information, 
then 5m b may be large at certain circumstances (e.g., Generally speaking, 9 2 is independent of energy E. Here, 



= 1-1 



r/2\ ( r 
^ exp l-2*r 



assuming a resonance with r—16 keV and -^- = 10 %, for simplicity, only s-wave case has been discussed. There- 
<5mb~59 %, 17 % at T = 0.1, 0.4 GK, respectively). Clearly, fore, it is sufficient to use the expression r a =ci\^EPi =0 (E)ocVEe^ 27TV , 
the variation 6mb exponentially depends on the stellar where the penetration factor Pi = o for the Coulomb bar- 
temperature, and therefore this assumption is not a good rier is approximated by a Gamow factor [3]. The quan- 
one for low temperatures. From the definition of Smb, an- tity c\ depends on the nuclear structure of the resonance, 
other complementary restriction for a narrow resonance The exponent is 2-kt)=ci\/ \ j E with C2—31 .2%Z\Zi^fJi de- 
should be expressed as r < j^jgg In 1 _g^f (r in unit of pending on the projectile-target combination. Note only in 
MeV). For instance, with 5mb=§ % and r/En<10 %, an- this equation the center-of-mass energy E is given in units 
other restriction for width is T < 0.017 Tg MeV. It should of keV. On the other hand, 7-ray partial widths depend 
be noted that in some literature a definition of 'narrow' only weakly on the energy via r i =c^E 2L+1 (x(E + Q) 2L+1 , 
was r<^En, but its ambiguity may give rise to certain where E 1 and L are the energy and multipolarity, respec- 



misleading in some sense. 



tively, of the 7-ray transition under consideration, and Q 
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is the reaction Q-value. The quantity C3 depends on the obviously this is not a good assumption for certain circum- 

nuclear structure of the two levels involved in the inter- stances (see discussion for the FIRST assumption). As a 

action. Therefore, for the usual particle-gamma reaction conclusion, except the SECOND reasonable assumption, 

channel (e.g., (p, 7) or (0,7)), the variation of quantity the FIRST and THIRD assumptions only hold for certain 

r a r y (oc^/E(E + Q) 2L+l e~ C2 ^ 1 / E ) over the resonance re- circumstances. 

gion can be expressed as Now, let's turn to Equ. IS Here an approximation (i.e., 

, / Af] \ 2L + 1 arc tan(-^ft)~-|) was utilized, and hence one can get 

xexp^U^W- 1 ). j o {E _ ER )2 + {r/2? dE -Yj2- 

Here, AE defines the energy interval between the reso- This approximation causes about a 3% error for a con- 
nance energy Er and the energies where the Maxwell- dition of r/E R — 10%, and thus it can be regarded as a 
Boltzmann function crossing the Breit-Wigner function, reasonable one. The integration ranges, from to 00, were 
It is obvious that the dominate variation takes place in used in Equ. [5J and therefore the quantity r a r b can't be 
the third exponential term. Given a condition of AE/Er considered as a constant and taken outside the integral. 
= 5 %, for p(p, 7) reaction, 5r a r 7 are only about 8%, The quantity r a r b can be roughly considered as a con- 
5%, 4% for resonances at Er=0.1, 0.5, 1.0 MeV, respec- stant over a very narrow energy interval, and therefore 
tively; while for 12 C(p, 7) reaction, 8r a r~, are respectively the integration ranges on Equ. O should be physically re- 
increased to 58%, 24%, 18% for the above resonance en- stricted [5]. A general equation has been derived here, 

ergies. Clearly the variation 5r a r^ is very sensitive to the „E a +n^ ^ 2 

/ (w- w \2, (r /9\2 dE = 7v^ arctan ( n )- ( 8 ) 

projectile-target combination via quantity C2- Therefore, JE R -n^ 1-° ^R) +l i / z J 1 / z 

the THIRD assumption is appropriate only for certain Ratiog between the values of quantity and those of 

projectile-target combinations for a resonance with J T/^<10 quantity _2_ a rctan(n) are 0.5, 0.7, 0.8 for n=l,2,4, respec- 

tively. This effect can be seen in the following Table [T] as 

Actually, the RRR per particle pair for a narrow reso- welL xhus for a small value of ^ the %ok reliable > deriva _ 

nance can be directly written as (use Eqs.l&E]) tion of Equ ^ is ma pp ropr i a te in the integration-range 

h 2 uj point of view; while for a large n, although the difference 



(av) = v2ir 



(nkTf/ 2 

3o / j? \ r r becomes small, the above-mentioned assumptions might 

x / exp(--g-) - " 2 b , 2 dE.(7) 

\ kl J [E — Er) + [1 J I) kg inappropriate in the 'narrow' point of view. Therefore, 







If one assumes quantity exp(—E/kT) as a constant over the classical analytic Eqs.rjJIHlare suggested not to be uti- 
the resonance region, Equ. [7] is equivalent to Equ. [6] But, lized in the future RRR calculations without caution. 
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3 Numerical Integration as, (i.e., Equ. (4.59)) in [3], 

a ( E )BW = a R 



b _e r r p {E) r,{E) (r R /2f 



e r p (E R ) r^En) (e - E R y + [r(E)/2]* ■ 

In the past, some researchers utilized the numerical in- (9) 

tegration of Equ. Q] with a broad-resonance formula to The energy dependence of the proton and 7-ray partial 

calculate the resonant-reaction-rate (RRR) of a narrow res- widths are given by 

onance, while others still utilized the simple analytic ex- ^\ rtr \ P(=o{E) 

pression of Equ.JSJ In order to check their differences, as an 

( E Qp~i — \ 

example, the RRR of a key stellar reaction 21 Na(p,7) 22 Mg ^(E) = r y (E R )B 7 y ^ R + Q — — E f ) ' 

is discussed below. There is a known resonant state in the where ig the (p; ^ react ion Q-value (Q=5.508 MeV [6]), 

compound nucleus 22 Mg at £ fl =0.821 MeV ( J*=l+, £=0) ^ is the primary 7 _ ray branching ratio to the final state 

with a proton width of T p =16 keV and a resonant strength at ^ (Her6j ^ ig assumed to be unity and Exf=Q for 
ofw 7 =0.556eV (here r 7 «^=1.48eV, and hence r«r p ) [6]. the ground gtate )_ The penetration factors, P t , are cal- 

The corresponding 7-transition to the ground state in cuiated by using a RCWFN g] code with a c h an nel-radius 

22 Mg is a pure All transition (i.e., 1^ >0^) with a mul- p a rameter of do = l 35 fm 

tipolarity L=l. According to the previous definition, this The numerically and ana l yt ically calculated resonant- 
resonance {r/E R ~2%) can be 'reasonably' treated as a rea ction-rates (RRRs) are compared in Table 1, where the 
narrow one. corresponding ratios were listed for different temperatures 
Figure 1 shows the amplitudes of relevant functions and integration ranges used in Equ. [7J It can be seen that 
used in the equations for this resonance state at the 'Gamow- two methods give almost the same results at 'Gamow- 
peak' temperature T=1.6 GK. The narrow- resonance cross peak' temperature T=1.6 GK for a full-range integration 
section a RW (E) is calculated by assuming the constant (0— >oo). For narrow integration ranges, the analytic re- 
partial widths r p , r y (and r). Obviously the tenden- suits are larger than the numerical ones. This shows that 
cies for this cross section and the corresponding integrand the integration should be computed for a full range, and 
curve (indicated by a grey thick solid line) are physically it is not a difficult task anymore with current computers, 
unreal towards lower energies. The reason is the assump- Notably, at lower temperatures, the numerical results be- 
tion of constant widths is invalid out of the resonance come much larger than the analytic ones (see the 2nd col- 
region. If a broad-resonance formula is used in the cal- umn in Table 1). This implies a higher-energy resonance 
culation, the integrand curve (indicated by a thick solid can contribute to the reaction rate at low temperature sig- 
line) become physically reasonable at low energies. Here, nificantly (or tremendously) more than the previous ana- 
the broad-resonance cross section cr RW {E) can be written lytical result. Furthermore, presume this -Er=0.821 MeV 
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resonance locating at 0.160, 0.323 MeV, the calculated Two years ago, a new book "Nuclear Physics of Stars" [7] 
Ratios for a full-range integration are 0.78, 0.90, respec- was published, where the definition of a narrow-resonance 
tively, at the relevant 'Gamow-peak' temperatures of 0.14, seems reasonable. But, we should realize that a question 
0.4 GK. This conclusion demonstrates that the previous is hidden in this sort of calculation, i.e., is it physically 
analytic results may be underestimated by a considerable appropriate to extend the broad-resonance Breit-Wigner 
amount for a resonance with a several keV width even formulism (esp. for a narrow resonance) down to very low 
at their 'Gamow-peak' temperatures. Although the con- energy? since the energy gap between the resonance and 
dition r/Eji<10 % still holds in these cases, Smb varies the Gamow-peak is quite large, it's more than hundreds 
from 6 % to 37 % for the temperatures listed in Table 1. or thousands of resonance width sometimes. However, we 
It indicates that, in the present case, the resonance width can't address this question in the present work. In ad- 
is not narrow enough to use the analytic equations at low dition, some authors [51|7j thought that the product of 
temperatures, and hence another restriction on the width Maxwell-Boltzmann distribution and cross section could 
should be set, e.g., r< 3.4 keV for a temperature of 0.2 give rise to another maximum caused by the low-energy 
GK, and it is about five times narrower than the experi- wing of a resonance but, our calculation shows the appear- 
mcntal value (_T=16 keV). ance of this low-energy maximum depends on the temper- 

ature condition, for instance, there is no such maximum 
4 Conclusions (but a plateau) in our example beyond T^0.5 GK (of 

course, no maximum at T=1.6 GK in Fig. [Tj). 



The validity of the resonant-reaction-rate (RRR) equations 
for an isolated narrow resonance have been reexamined, 
and it reminds us not to use those analytic equations with- 
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